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Abstract. Given any finite set T of (n — 1)- dimensional subspaces of R" we give 
examples of nongaussian probability measures in R" whose marginal distribu- 
tion in each subspace from J- is gaussian. However, if J 7 is an infinite family of 
such (n — 1) -dimensional subspaces then such a nongaussian probability mea- 
sure in R" does not exist. 
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1 Introduction 

Starting with the simple example of E. Nelson as cited by W. Feller in JTJ 
we have from the papers of B.K. Kale |3|, G.G. Hamedani and M.N. Tata 
|2| and Y. Shao and M. Zhou 01 etc., as well as Section 10 of J. Stoyanov's 
book (51, several examples of bivariate and multivariate nongaussian dis- 
tributions under which many linear functionals can have a gaussian dis- 
tribution on the real line. These results suggest the possibility of char- 
acterizing a gaussian distribution in R" through properties of classes of 
linear functionals. Motivated by Nelson's example in [1 j and the bivariate 
construction in we introduce a perturbation of the standard gaussian 
density function in R" exhibiting the following interesting features: (1) 
Given any finite set {Sj, 1 < j < N} of (n — 1) -dimensional subspaces 
it has a marginal density function which is standard gaussian in each Sj, 
j G {1,2, ...,N}; (2) There can exist linear functionals whose distributions 
may have nonunimodal density functions; (3) For certain choices of sub- 
spaces the nongaussian perturbation can be so chosen that any real sym- 
metric measurable function of all the n coordinates has its distribution pre- 
served. In particular, the sum of squares of all the coordinates can have the 
X 2 distribution with n degrees of freedom. 



We also demonstrate the following characterization of the multivariate 
gaussian distribution. Suppose {Sj,j = 1,2,...} is a countably infinite set 
of (n — 1) -dimensional subspaces of R" and u is a probability measure in 
R" such that the projection of w in each subspace Sy is gaussian. Then ft 
itself is gaussian. This is a generalization of the characterization in [2] and 
a more precise version of the result in [4j. 

Our proofs follow the steps in [2] and use some additional geometric 
and topological arguments of a very elementary kind. 

2 A perturbation of the gaussian 
characteristic function 

We begin by examining a small perturbation of the characteristic function 
of the ft-variate standard gaussian distribution with mean vector and 
covariance matrix I as follows. Choose and fix any homogeneous polyno- 
mial V of even degree 2k in n real variables t\, t^, t n and define 

<$>{p,E,a,V) = e-iW +e e-$°*Wv(t), t e R" (2.1) 

where t = (t\, t n ) T , £ is a real parameter and a is a parameter satisfying 
< a < 1. Here 

\t\ 2 =(tl + ... + t 2 n ). 
Clearly 0(-; e, a, V) is a real analytic function on R" satisfying 



O(0;e,(r,P) 
<P{-t;e,a,r) 

Let 

Z v = {t\V(t) 



= X 

= ®(t;e,<r,V). (2.2) 
= 0,fGR"} (2.3) 



be the set of zeros of V in R" . 

Since we are interested in the inverse Fourier transform of <3> we intro- 
duce the renormalized polynomial : V : in the form of a formal definition. 

Definition 2.1. Let 



and let H m (x) be the m-th Hermite polynomial defined by 
d m 

— m{x) = {-l) m H m {x)Vl(x), m = 0,1,2,... 
(as in Feller /[!][). For any real polynomial V in n real variables given by 

T>(+« ir, i \ — n f m U™2 j.m„ 

I l Li l n) — u m\,mi,...,m n l \ t 2 " n 

m 

its renormalized version :V:is defined by 

■ T '• {X\,:-)Xfi) = ^mi,m2,...,m n timi{^l)H m2 ^X2)---Hm„{Xn)- 

m 

Note that for a homogeneous polynomial, its renormalized version need not be 
homogeneous. 

Since the function O in ( |2.1| is in Li (R n ) its inverse Fourier transform 
/ is defined by 

f(x;e,<r,V) = J ' e- itTx <S>{t;e,o-,V)dt 1 dt 2 ...dt n 

= _1 e 4W 2 + e ^ [ e-' tTx e- 1 ^ 2 V(t)dt 1 ...dt n . 

(2.4) 

First, we note that 

(27T)« y cr" V ^ / 

Repeated differentiation with respect to X\, x%, x n shows that for the ho- 
mogeneous polynomial V of degree 2k we have 

^ J e- ifx e-^\ t \ 2 V{t)dhdh...dt n 

£-(?)} 



C n V 9xi' ' dx n 



:V: ' e-^ 

cr' 1 ™ \a a J (v27r) n 



3 



Thus the inverse Fourier transform ( |2.4| | assumes the form 
f(x;e,(T,V) 

-li*i 2 (i + izigi , v , (* *l) r^PM)}. 

{V2tt)" \ <r n+lk \a' ' a J J 

(2.5) 

Since, by assumption, 1 — a 1 > the positive constant K(c, T 7 ) defined 



by 

K[*,V) = sup l^ : (^-^)l e -iW'(i-^) (2.6) 

is finite and for all x G R n 

f(x;e,a,V) > if |e| < K^ 1 (a,T') 
we observe that 0(-;e, c, T 5 ) is a real characteristic function of the proba- 



bility density function f(-;e,cr,V) defined by (2.5 for any 

e G [— K^ 1 (cr,V) , K^ 1 (cr,V)}. Here we have made use of property (2.2). 

Thus we can summarize the discussion above as a theorem. 

Theorem 2.2. Let < a < 1, V be a real homogeneous polynomial in n 
variables of even degree 2k, K(a,V) the positive constant defined by \2.6\ a nd 
ee [-K" 1 (cr, V), K _1 (cr, V)]. Then the function Q(-; e, cr,V) defined by \2.1} is 
the c haracteristic function of a probability density function f(-;e,a,V) defined by 



(2.51. Under this density function f(-;e,a,V) the linear functional x i — > a T x 
with \a\ = 1 has characteristic function (p a and probability density function g a 
on the real line given respectively by 

<p a (t) = e-i t2 +e V(a)e- 1 2 a2f2 t 2k , t € R (2.7) 



wW _-^{^ + H^W^(^).-^}. (2 .s, 

In particular, for any a G Z-p, the linear functional a T x has the normal 
distribution with mean and variance \a\ 2 but f(-;e,a,V) is a nongaussian 
density function for any t E [— K~ x (a, V),K~ 1 (ct, V)] \ {0}. 
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Proof. The first part is immediate from the discussion preceding the state- 
ment of the theorem. To prove the second part we note that the character- 
istic function cp a (t) of the linear functional a T x under the density function 
f(-;e,cr,V) is e,cr,V) and (2.7 follows from (2.1 1 and the homogene- 
ity of V. Now (|2.8[l follows from Fourier inversion of (|2.7|. If 7^ a G 



thenO = = P 



and therefore 



<P a{t) = e i* . 

W 

Hence a T x is normally distributed with mean and variance | a | 2 . □ 

Corollary 2.3. Let {Sj,l < j ' < N} be any finite set of (n — l)-dimensional 
subs-paces ofR n . Then there exists a nongaussian analytic probability density 
function whose projection on Sj is gaussianfor each j E {1, 2, N}. 

Proof. By adding one more (n — 1) -dimensional subspace to the collection 
{Sj, 1 < j < N}, if necessary we may assume without loss of generality 
that N is even. Choose a unit vector «W G Sj- for each j and define the 
homogeneous real polynomial V of degree N by 

N 

V{t) =Y\a {j)T t, t G R' ! . 

Clearly 

N 

V(t) = if f G [jSj. 

;'=i 

In other words 

N 

U s ; C Zp- 

If we choose ^ to be the probability measure with the density function 
f{-;e,a,V), ^ e in [-K- 1 (c7 / 7 7 ) / K- 1 ((7 / P)] in Theorem 2.2 it follows 
immediately from the last part of the theorem that every linear functional 
of the form b T x has a normal distribution with mean and variance \b\ 2 
whenever b G Zp. This completes the proof. □ 
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Remark 2.4. In the context of understanding the modes of the density function 
ga (x) given by (|2.8|> it is of interest to note that 



Indeed, this is obtained by straightforward differentiation and using the recur- 
rence relation H 2 k+i{x) = xH 2 k(x) — H' 2k (x). 

Example 2.5. Let n be even, 



V{h,t 2 ,...,t n ) 



ht 2 ...t„Yl(t^ - tj] 



= t\t 2 ...t n 



1 


1 


1 


A 






A 


A ■ 


■ A 


,2(n-l) 


2(n-l) 
l 2 


2(n-l) 



(2.9) 



Then V is a polynomial of even degree n 1 , which is antisymmetric in the 



variables t\, t 2 , —, t n . The renormalized version : V 



V : (x 1 ,x 2 ,...,x„) 



H 3 (*i) 



Hi(x 2 ) 
H 3 (x 2 ) 



H 2 n+l(*l) H 2 „ + i(x 2 ) 



of is given by 



Hi(x, ; ) 
H 3 (x n ) 



H 2n+ i(x n ] 



(2.10) 



In particular, : V : is antisymmetric in the variables Xi, x 2 , x n . Fixing 
< a < 1 we get for each e G [-K.- l {a,V),K.- 1 {a l V)] l with K{a,V) being 
determined by (2.6 1, (2.10' and /c = \n 2 , the probability density function 
f(-;e,cr,V) given by 
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f{x;e,a,V) 



(y/2n)> 



r »(n+l) 



B 



2n+l 



Him. 



H 



2n+l 



■X2' 

- a ■ 



Him, 
Him, 



H 



2n+l 



m 



(2.11) 



By Theorem 2.2 and its Corollary we conclude that the projection of 
this density function on the (n — 1) -dimensional hyperplanes {x\xj = 0}, 



1 < j < n; {x\x( 



0}, 1 < i < j < n and {x\x t + x ] = 0}, 1 < i < j < 



n are all (n — 1) -dimensional gaussian densities. 

If g{x\,Xi, ...,Xn) is any bounded continuous function which is sym- 
metric in the variables X\, x 2 , x n then the function g : V : is an antisym- 
metric function in R" and therefore 



f 1 |y|2 

/ (g:V :){xi,X2,...,x n )e ^ 1 dx 1 dx 2 ...dx n = 0. 

JR" 

Thus 

/ g(x lr x n )f(x; e, a, V)dx 1 dx 2 ...dx n 

JR" 

f 1 1 2 

= / g(x 1 ,x 2 ,...,x n ) — e^ x \ 2 dxidx 2 ...dx n . 

In other words, for 0^£G [— K _1 (a,V), K^ 1 (a, V)], any symmetric mea- 
surable function g on R" has the property that its distribution under the 
nongaussian density function f(x;e,a,V) in (2.11 1 is the same as its dis- 
tribution under the standard gaussian density function with mean and 
covariance matrix I. 

Example 2.6. We now specialize Example 2.5 to the case n = 2, a = 2~ 1,/2 
when 

V{h,t 2 ) = hh(t\-t\), 
:V:(x lt x 2 ) = H 3 (x 1 )H 1 (x 2 )-H 1 (x 1 )H 3 (x 2 ) 
= x\x 2 — x\x\. 
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A simple computation shows that 



K(cr,V) = 8sup|xfx 2 -x|xi| e~^ +x ^ 
= 128 e~ 2 . 

This supremum is easily evaluated by switching over to the polar coordi- 
nates x\ = rcosO, X2 = rsinO. Then 

f(x;e,a,V) = -U"^* 2 ) {l +32e(^ 2 - x 3 2 xi)e-^ x * +x %} (2.12) 
which is a probability density function whenever 

£ < • 

11 - 128 

At £ = 0, it is the standard normal density function with mean and 
covariance matrix I. We write rj = 32 £ and express the density function 
d2l2|as 



f, 1 {x l ,x 1 ) = JL g -i(*?+*i) {l + j/^f^-x^Oe-H^+^J (2.13) 
where (See Fig. Q.) 

When a = (sin9,cos9) the density function go of the linear functional 
x i — > X\sin9 + X2COS0, under is given by the formula ( |2.8[ l of Theorem 
2.2 as 

= |l " ^1^M(4^ _ I2x 2 + 3) e -^ 2 | . (2.14) 



Thus 



^(x) = —=e~^ x { ez* - 1 — '—— — 2 — i.(4 x 4 _ 20x 2 + 15)e 2* 



It is not difficult to find values of rj in the range (0, \e 2 \ and angle for 
which 



el* 2 - v ll^ n M(4x 4 - 20x 2 + 15) = 0,x ^ j ^ 0. (2.15) 

This reveals the possibility of nonunimodality of the density of some linear 
functionals under the joint density fa. For an illutration cf. Fig. |2]). 



8 



Figure 1: Bivariate density f n (xi,X2) at r\ = e 2 /4. 

3 A characterization of gaussian distributions in R n 

In the context of the Corollary to Theorem 2.2 we have the following char- 
acterization of a gaussian distribution in R" when the number N of (n — 
1) -dimensional subspaces in the corollary is countably infinite. 

Theorem 3.1. Let {Sj,j = 1, 2...} be a countably infinite set of (n — 1)- dimen- 
sional subspaces o/R" and let jibe a probability measure in K" whose projection 
on Sj is gaussian for each j = 1,2,... . Then ]i is gaussian. 

Proof The fact that the projection of ji on the two distinct (n — 1 )-dimensional 
subspaces Si and S2 are gaussian implies that the multivariate Laplace 
transform y. of ]i given by 

ji(zi,...,z n ) = / exp(ziX! + ... -\- z n x„)u(dx 1 dx 2 ...dx n ) (3.1) 
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9 = 0.57° 



9 = 6.87° 





9 = 14.32° 



9 = 54.43° 





Figure 2: Nonunimodality of gg. 

is well-defined for z S C" and analytic in each of the complex variables Zj, 
j = \,...,n. Let m and E be respectively the mean vector and covariance 
matrix of the K" valued random variable x with distribution ji. 

Choose and fix a unit vector e Sy- for each = 1, 2, ... . Suppose 



,0') r = 



Since 



(«;'!/ •••/«;«)/ / — 1/2, 



a,- = max a,> . 

y l<r<« ; 



E4 



1, V; 



r=l 



it follows that Kj > n 1/2 , V There exists an tq such that fl ; > = ay for 
infinitely many values of ;'. Restricting ourselves to this infinite set of fs 
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and assuming tq = 1 without loss of generality we may as well assume 
that 



— ifljli •••/ "in) / 
/i I = max |fl,>| V /' = 1,2,..., 

7 KKn 7 



J /l| 



> n 



-1/2 



V;. 



Now consider the (n — 1)- dimensional vector b^ ; defined by 



«/2 A;'3 



4 /» 



flyi fljX "/I 



, 7 = 1,2,. 



where 



*,1 



< 1 Vr = 2,3, ...,n. 



Thus all the vectors b^ are distinct and they constitute a bounded count- 
able set in R*" -1 ). Define the set 

D = P| |s|s G R^ 1 ', (b^ - b (0 ) T s ^ 0} . 

/<i 

Being a countable intersection of dense open sets it follows from the Baire 
category theorem that D is dense in R("~ 1; . Let now 

s = (s 2 ,s 3 ,...,s n ) T 

be any fixed point in D. Define 

s/i = -b^ T s,j = 1,2,.... 

By the definition of D, {syi, / = 1, 2, ... } is a bounded and countably infinite 
set of distinct points on the real line. Furthermore 

0/lS/l + A/2S2 + - + a jnSn = V /'. 

In other words, (sji,S2, —,s n ) T G Sj for each By hypothesis the linear 
functional has a normal distribution with mean 

Sjittii +S2?«2 + ... + s„m n and variance (syi,S2, ...,s„)Z(syi,S2, s„) T . Defin- 
ing 

j/>(zi, ...,z n ) = exp(m T z + -z T Zz), z G C" 
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we conclude that the Laplace transform ji defined by (3.1 1 and the function 
xp satisfy the relation 

fi(Sj 1 ,S2,...,S n ) = tp{Sj 1 ,S 2 ,...,S„) 

for j = 1,2,... . Since fi(z,s 2 , ...,s n ) and ip(z, S2,...,s„) are analytic functions 
of z in the whole complex plane and they agree on the infinite bounded 
set {sj\,j = 1,2, ...} it follows that 

g(z,s 2 ,...,s„) = f(z,s 2 ,...,s„)Vz G C. 

Since this holds for all (S2, ...,s„) T £ D which is dense in R^" -1 ) and both 
sides of the equation are continuous on R" we have 

f{si,s 2 ,...,s n ) = ip(si,s 2 ,...,s n ) 

for all (si, s 2 , s„) T € R". This implies that \i is a gaussian measure with 
mean vector m and covariance matrix Z. □ 
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